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Problem  StaLemciit 


f 


1.  Tnirodii.' t  ion 


J 


The  Air  Force  lurrently  stocks  and  manages  i  total  of  1.48 
million  expendable  line  items  in  its  multi-level  inventory  system. 

The  dollar  investment  in  inventories  of  these  relatively  low  cost 
items  is  over  1.1  billion.^ 

Efforts  have  been  directed  toward  developing  and  applying  models 
as  approximations  to  reality  in  order  to  manage  these  large  inven- 
tories.  These  models  allow  managers  to  conceptualize  the  nature  of 
the  supply  environment,  the  inventory  system,  the  various  subsystems 
within  tlie  inventory  system,  and  the  relationships  between  environment, 
system  and  subsystems.  An  important  benefit  of  this  modeling  approach 
is  the  fact  tliat  an  appropriate  model  can  provide  a  moans  of  minimizing 
total  inventory  costs,  and  yet  allow  the  system  to  remain  responsive  to 
the  needs  of  the  item  users. 

The  Air  Force  curreYitly  uses  a  modified  version  of  the  classic 
Wilson  lot  size  economic  order  quantity  (EOQ^  model  for  inventory 
system  management  of  expendable  line  items.  The  objective  of  this 
model,  designated  the  D062  model,  i.s  to  minimize  total  variable  costs 
of  ordering  items  and  maintaining  them  in  inventories.  ;  The  major 
constraints  are  budget  restrictions  and  the  requirement  for  adequate 

•  r 

response  to  the  needs  of  the  item  users. 
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I  ' 

Several  s  iy.n  i  f  leant  Factors  of  the  D()62  moilel  liave  not  been  n^le- 

quately  refined  and  have  caused  problems  in  the  practical  application 

of  the  model  to  inventories  within  the  Air  Fore?  supply  environment. 

Major  problems  include  inadequate  estimates  of  neaningful  cost  to 

order  and  cost  to  liold  factors,  problems  of  bud>ot  constraints,  and 

•> 

lnae<-(irate  forecasts  of  future  item  demand  patt?rns.“ 

Some  of  those  problems,  and  others  not  stated,  were  addressed  in 
the  development  of  the  ne\^;  A022  liOCJ  model,  whicn  the  Air  Farce  intends 
to  adopt  as  a  replacement  to  the  D062.  Many  research  efforts  have 
been  initiated  to  make  the  practical  application  of  this  nov  model 
more- refined  and  therefore  more  effective.  This  paper  is  the  result 
of  one  such  research  effort.  The  area  of  concern  addressed  by  this 
research  effort  was  the  problem  of  inaccurate  forecasts  of  future 
item  demand  patterns. 

Background 

The  need  for  routine  forecasting  methods  exists  at  all  levels  in 
the  Air  Force  multi-level  inventory  system.  Air  Force  base  level 
forecasting  is  particularly  important  since  requirements  at  this  level 
drive  total  system  stock  levels.  There  are  two  primary  objectives 
at  base  level  for  the  EOQ  management  model.  The  first  is  to  prevent 
stockouts  and  subsequent;  frustration  of  user  needs.  The  second  is  to 
manage  base  stock  levels  to  minimize  the  sum  of  the  inventory  carrying 
costs  and  ordering  costs.  Accurate  forecasts  of  base  level  item  demand 
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patterns  as  an  input  to  tlio  l^OQ  model  improves  tie  model's  capability 
to  meet  these  objectives.  In  turn,  base  level  demand  patterns  are 

inputs  into  hip, her  echelon  niOQ  management  systenu,  which  also  must 

! 

resolve  probli.'ms  of  stockouts  versus  inventory  costs.  The  key  Chen, 
to  improving  totkl  system  response  in  tlie  area  of  stockouts  versus 
inventory  costs  is  more  accurate  base  level  forecasting  of  item 
demand.  ‘ 

Purpose  of  Research  ^ 

\ 

The  present  method  of  unweighted  moving, averages  to  forecast  base 
level  demand  patterns  has  a  tendency  to  over-stock  inventories  during 
periods  of  declining  demand  and  under-stock  inventories  during  periods 
of  significant  demand  increases.  This  over-stocking  or  under-stocking 
of  inventories,  applied  to  Che  world  wide  scale  of  the  Air  Force  supply 
environment,  resul-ts  in  large,  unnecessary  inventory  holding  costs  on 
the  one  hand  and  emergency  procurement  (ordering)  costs  on  the  other. ^ 
The  purpose  of  this  research  effort  was  to  examine  and  compare 
five  methods  of  base  level  demand  forecasting  which  can  be  considered 
as  alternatives  to  the  present  method  of  unweighted  moving  averages. 

All  six  methods,  including  unweighted  moving  averages,  were  evaluated 
according  to  relative  accuracy  in  describing  the  actual  base  level 
demand  patterns  of  a  statistically  chosen  stratified  sample  of  items 
managed  at  base  level  under  the  Air  Force  EOQ  model. 
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This  study  v;as  exclusively  coucc-rned  with  e<pendnble  itcns  which 
are  stocked  based  upon  D062  EOO  procedures.  Expindable  itens  arc  those 
which  are  “consumed  in  use  or  which  lose  their  original  identity  during 
incorporation  ipto,  or  attachment  upon,  another  issenbly" . ^  These  items 
coranonly  referred  to  as  EOQ  items,  arc  designated  by  expcndabil i ty ,  ro- 
pairahility,  recoverability  category  (ERRC)  codes  ''X32"  and  "XB3" ,  among 
others.  ; 

This  study  was  strictly  limited  to  a  sample  of  XB2  and  XB3  items 
stocked  by  the  Air  Force,  and  present  in  a  substantial  number  of  base 
level  inventories.  '■■■The'  complete  definitions  of  XB2  and  XB3  ERJIC  coded 

items  are  as  follows:  ! 

.1 

XB2 :  Expendable,  nonrecoverable  (no  repair)  items  witli  a  projected 

annual  requirement  of  $10,000  or  more  regardless  of  unit  price. 

XB3;  Expendable,  nonrecoverable  (no  repair)  items  with  a  projected 
annual  requirement  of  $10,000  or  less  regardless  of  unit  price. ^ 

Five  basic  forecasting  methods  were  evaluated.  Tliese  methods  were: 
unweighted  moving  averages,  least  squares,  and  the  methods  of  single, 
double  and  triple  exponential  smoothing.  A  sixth  method  was  generated 
by  modifying  the  single  exponential  smoothing  method.  The  modification 
consisted  of  incorporating  an  adaptive  exponential  smoothing  constant 
which  can  be  changed  according  to  the  affects  of  exogenous  variables. 

The  results  obtained  with  this  sixth  method  were  also  compared  against 
the  results  of  the  five  basic  methods. 


The  validity  of  n  foroc.i  .t  i I'.y  -.i  t'.ioJ  ri-sts  wit  i  iis  ability  to 
perform  in  a  real  world  oiivi  ronnent .  I'iierefore  it  was  -iecided  to  test 
Che  six  given  forecasting  methods  u..ng  actual  <.  emand  experience.  The 
objective,  then,  of  this  study  was  to  evaluate  the  given  forecasting 
mctliods  usiii;.’,  an, actual  sample  of  ecanonic  ordci  quantity  items  stocked 
at  a  base  level  consolidated  supply  activity,  llie  metiiods  were  com¬ 
pared  to  determine  if  significant  differences  existed  in  the  accuracies 
of  their  rcspet:tive  forecasts. 

Brief  Summary  of  Find ings 

This  summary  describes  the  tluoe  most  accurate  of  the  six  evaluated 
forecasting  methods  over  the  total  sample  set  of  316  items  and  throe 
subsets  of  tlio  total  set.  The  subs>!ts  are:  Subset  1,  the  71  items  with 
relatively  liigli  level,  non-erratic  demand  processes;  Subset  2,  the  199 
items  v/ith  relatively  low  level,  non-crratic  demand  processes;  and 
finally.  Subset  3,  the  '16  items  with  relatively  erratic  demand  processes. 
An  erratic  demand  process  means,  in  this  study,  that  no  discernible 
constant  level,  linear  trend  or  curved  trend  in  demand  v^as  present 
over  the  entire  15  month  period  for  a  given  item. 


Forecasting  Method 

Me.an  Squared 
Error 

Mean  Forecast 
Error 

Variance 

1. 

Single  Exponential 

Smoothing 

16.02 

1.74 

13.02 

2. 

Moving  Average 

17.84 

1.78 

14.72 

3. 

Double  Exponential 

Smooth ing 

23.89 

1.89 

20.39 

Table  1 .  Total  Sot  (316  Items) 


Forecast ing 

Method 

Mean  Squared 
_ Error 

Mean  Forecast 
Error 

Var iance 

1. 

Single 

Kxponent ial 

Smooth ing 

67.99 

5.61 

38.16 

2. 

Moving 

Average 

76.18 

5.68 

44. -W 

3. 

Double 

Exponent ial 

Smoothing 

103.19 

6.1A 

66.33 

Tnhle  2.  SuliseC  1  (71  Itens) 


For 

Mean 

oc. 1st  ing.  Method 

Squared 

Error 

Mean  Forecast 
Error 

Var iance 

1. 

Single  Exponential  Smoothing 

.36 

.49 

.12 

2  ^ 

Moving  Average 

.37 

.49 

.12 

3. 

Double  Exponential  Smoothing 

.40 

.50 

.14 

Table  3.  Subset  2 

(199  Items) 

Mean 

Forecasting  Method 

Squared 

Error 

Mean  Forecast 
Error 

Var  iance 

1. 

Double  Exponential  Smoothing 

3.56 

1.32 

1.36 

2. 

Moving  Average 

3.63 

1.31 

2.01 

3. 

Single  Exponential  Smoothing 

3.81 

1.35 

2.04 

Table  4.  Subset  3 

(46  Items) 

Single  e;cponential  smoothing  provided  the  most  accurate  forecasts 
for  the  total  set  of  316  items  and  also  for  component  subsets  1  and  2,  based 
on  rankings  according  to  variance,  mean  forecast  error,  and  mean  squared 


error.  Only  oa  Lhe  items  with  relatively  erratic  tleraand  processes  did 

I 

the  order  of  aemracy  Chanp.e.  In  this  case,  doiible  exponential  smoothing 

was  the  most  accurate  method.  The  moving  averai;e  method  consistently 

I 

ranked  second  for  tlie  total  set  and  subsets  1,  ?.  and  3. 

However,  a-fcording  to  the  Kruskal-Mal  1  is  t.'St  which  coriipared  the 
absolute  average  error  distributions,  no  signifreant  difference  was 
detected  in  tlie  relative  effectiveness  of  the  six  forecasting  n>ethods.  \ 
confidence  interval  of  95  percent  was  used  for  this  test. 

Organization  f 

The  remaining  ttirce  chapters  of  this  research  paper  are  concerned 
with  the  methodology  used  in  the  study,  an  analysis  of  results,  ,and  a 
detailed  summary.  Chapter  II.  Methodology ,  contains  a  description  of 
the  data,  a  presentation  of  the  forecasting  methods  and  a  discussion  of 
the  experimental  design.  Chapter  *III.  Analysis  of  lesiilts,  compares 
forecast  method  performance.  Chapter  IV.  Detailed  Summary,  presents  a 
detailed  summary  of  the  findings,  lists  the  limitations  inherent  in  the 
research  methodology,  and  provides  recommendations  for  action. 
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Descr ipt ion  of  Data 

All  six  methods  of  forecasting  were  evaluated  according  to  rela¬ 
tive  rofincment  and  accuracy  in  describing  the  actual  deinand  patterns 
of  a  stat  i.st  ically  chosen  stratified  sample  of  items  managed  at  base 
level  under  Che  Air  Force  KOO  model.  The  methodology  used  in  clioosing 

the  sample  is  important  because  the  results  generated  by- testing  the 

» 

six  forecasting  methods  apply  only  to  the  demand  processes  of  wliich 
the  sample  was  representative.  This  section  is  therefore  concerned 
willi  the  source  of  data  in  the  sample,  the  techniques  used  In  sample 
selection  and  a  description  of  how  the  sample  was  used  in  the  actual 
testing  of  the  six  forecasting  metliods. 

The  only  source  of  detailed  base  level  demand  data  for  bOQ  man¬ 
aged  items  is  the  consolidated  and  daily  transaction  registers  main¬ 
tained  at  base  consolidated  supply  activities.  Such  registers  are 
detailed  records  of  all  transactions  on  all  active  i terns'  stocked  at 
the  base.  Registers  are  maintained  for  the  previous  and  current  cal¬ 
endar  years  so  that  historical  data  is  available  for  one  to  two  years 
back  depending  on  the  current  calendar  month. ^ 

Source  of  Data.  The  empirical  basis  for  this  research  effort 
was  15  months  of  demand  data  (I  January  1973  through  31  March  1974) 
for  selected  expendable  supply  items.  The  data  source  for  the  selec¬ 
ted  items  was  the  transaction  registers  maintained  by  the  21  Supply 
Squadron,  Elmendorf  Air  Force  Base,  Alaska. 
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Time  Periotls  for  Dem.iiul  Djl.i.  Tlie  actual  denand  data  i.-as  aj'.j'.re- 
gated  into  15  monthly  tota^s  for  eacli  item.  The  first  10  months  of 

,  I 

demand  Iiistory  (1  Janiiarjf  1973  through  31  October  1973)  were  used  as 
the  experience  base  from  which  to  make  monthly  forecasts  for  the  fol¬ 
lowing  5  months  (i  November  1973  through  31  March  1974).'  Th.e  10 
month  experience  base  was  designated  as  the  base  period.;  The  follow¬ 
ing  5  months,  then,  was  designated  the  forecast  period.  -  The  purpose 
of  the  actual  demand  data  in  the  forecast  period  was  to  use  it  for 

error  comparison  in  evaluating  the  relative  accuracy  of  <the  given 

i 

six  forecasting  methods,  i  * 

Sampling  Technique.  The  initial  sample  included  al,l  XB2  and 
XB3  items  for  three  Air  ^orce  weapons/support  systems.  Those  items 
which  had  at  least  one  demand  during  the  two  year  period  1  July  1972 
through  31  June  1974  were  included  in  the  sample.  The  selection 
process  for  the  three  systems  was  essentially  random,  wi,th  the 
criteria  that  the  systems  be  commonly  deployed  at  many  bases,  and 
not  unique  to  the  Alaskan  theater.  ; 

The  sample  items  for  the  three  chosen  systems  were  stratified 

I 

according  the  "ABC"  technique.  System  identity  was  not  maintained 

in  this  stratification  process.  The  final  sample  (selected  from 

1 

the  initial  sample)  was  chosen  to  be  reasonably  representative  of 
the  stratified  "B"  cost  group  of  all  XB2  and  XB3  expendable  items 
used  in  support  of  common  Air  Force  weapons/support  systems.  A 
brief  explanation  of  the  ABC  technique  is  pertinent  at  this  point. 


V 
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The  ABC  Toc’.hni<int;.  Inventory  manaf.amenc  invjlves  the  n.inageiiont 
of  individual  items.  Unless  each  item  stocked  is  under  a  reasonable 
degree  of  control,  the  aggregate  will  not  be  under  adequate  control, 

A  technique  is  needed  that  will  isolate  those  itens  that  have  a  high 

» 

total  dollar  inventory  turnover  per  year  from  these  that  have  a  rela¬ 
tively  low  total  dollar  ttirnover. 

The  high  dollar  turnover  items  should  be  of  primary; concern  to 
the  inventory  management  system  and  should  be  controlled  with  more 
precision  compared  to  low  dollar  turnover  items.  This  is  simply  an 
application  of  Pareto's  Principle  of  Maldistribution,  which  has  been 
expressed  as  follows;  "Very  often  a  small  number  of  important  items 
dominate  the  results  while  at  the  other  end  of  the  line  are  a  large 
number  of  items  whose  volume  is  so  small  that  they  have  little  effect 
on  the  results".^ 

Normally,  input  data  for  each  item  undergoing  ABC  analysis  should 
include  the  following: 

(1)  Stock  Number 

(2)  Unit  Cost 

(3)  Usage  per  time  period 

These  data  are  manipulated  in  the  following  manner: 

(1)  Multiply  usage  per  time  period  by  unit  cost  to  obtain  the 

I 

dollar  demand  for  each  item. 

(2)  Sort  the  items  according  to  dollar  demand  in  a  descending 
sequence. 
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(  i)  l-ist  tills  spquo.nct!.  Obt.'i'in  a  ctir:i:lat  iv  i’  liollar  JfT-.and  fii;- 

O 

lire  as  each  item  is  added  to  the  list.  , 

I 

The  listing  will  normally  showi  that  a  relatively  f<Jw  icens  have 

a  high  impact  on  us.ige  value.  Tlie  "A"  group  is  the  classif  icat  ion  of 

i 

highest  ilollar  demand  items.  It  is  usually  chosen  ns  tliose  items 
that  make  up  35  to  83  percent  of  the  total  cumulative  dollar  demand 
figure,  going  dov/n  the  list  from  the  highest  dollar  demand  item 
toward  the  lowest.  Usually  the  A  group  includes  only  aoout  20  per¬ 
cent  of  the  total  number  of  items. t 

I 

The  "B"  group  classification  is  usually  chosen  to  lie  95  to  98 

1 

percent  of  the  total  cumulative  dollar  demand  figure,  the  B  group 

i 

includes  those  items  starting  at  tiie  top  of  the  list  and  going  dowTi 
until  95  to  98  percent  of  the  cumulative  dollar  demand  is  reached. 

'file  B  group  adds  another  20  percent  of  the  total  items  to  those  selec¬ 
ted  by  the  A  group,  for  a  total  of  about  40  percent  of  clic  total  items. 
The  remaining,  low  usage  items  are  classified  as  the  "C.-  group,  which 
contains  about  60  percent  of  the  total  items. 

The  ABC  Tcclinigue  Applied.  The  initial  sample  of  items  chosen 
for  this  analysis  included,  as  previously  stated,  all  XB2  and  Xii3 
items  required  for  three  weapons/support  systems  in  common  Air  Force 
use.  For  inclusion  in  the  initial  sample,  each  item  must  have  had  at 
least  one  demand  during  the  two  year  period  1  July  1972  through 
31  June  1974.  The  three  randomly  chosen  weapons/support  systems  were: 
the  F4E  aircraft,  the  C-130  aircraft  and  the  Digital  Subscriber  Termi¬ 
nal  Kcpii lament  (DSTF). 
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Fol  low  in;.’,  till’  AlU:  1 1  i  hu  iiii.e,  .the  t  i’r,u  irpd  ;:-pit  liai.i  for  each  iLcm 
of  tho  initial  san'.'plc  were  rea-.n  pii  I  at  cd  in  tfib  fc  llowine  manner: 

(1)  The  dollar  demands  oL  the  initial  samide  of  XT>2  and  .\33  items 

were  calculated  .iccordinj'  to:  i 

=  L''j]  lar  demand  for  item  i  (i  =  1  , 2 , 3 , .  .  ,n) . 

Dj  =  Vj  X  P^,  where 

»•  I 

Vi  =  demand  volume, 

year 

P,  =  unit  price  and  ' 

n  =  total  number  of  items. 

(2)  All  items  v/ere  then  sorted  from  hijdiest  to  lowi’st  value  of  Di* 

(3)  After  listinj;  the  sorted  items,  the  Di  for  e.ich  item  was  .iddeti, 
in  sequence,  cumulatively. 

n 

^  D.  =  Total  cumulative  dollar  demand  per  vear. 

i^  " 

(A)  Items  selected  for  the  final  sample  included  those  starting 
with  the  highest  value  of  Di,  raovin;’  continuously  down  the  list  until 
95  percent  of  the  total  cumulative  dollar  demand  figure  was  reached. 

The  final  sample  was  therefore  the  B  cost  group  of  the  initial 
sample.  This  application  of  the  ABC  technique  selected  the  most  active 
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X!52  ami  XtVj  i  t  cn!:: ,  ilolla?'  v.isc,  in  i ‘ii'  i))  •"i-miar !  Ain  invoacorv  ;\'r  ;  iie 

.  •  1 

tlircc  rnndoinly  clin.son  sys  t.  duiin;’.  the  period  1  July  1472  il'.r.m.'.li 
31  June  197/i.  demand  tkrta  for  the  fir.al  sample  was  avail¬ 

able  for  the  period  1  January  1973  tliroupji  31  March  1974.  These  demand 
data  for  tlte  final  sample  of  items  were  used  to  evaluate  ti\e  six  fore¬ 
casting  methods  considered  in  this  researcli  effart.  Thcro  were  977 
items  in  the  initi.’.l  .sample.  The  B  <;roup  redticed  the  final  sample 
size  to  316  i  lotus,  approximately  32  percent  of  the  total. 


1^1  r CO astinp.  Me f  hod s 

The  objective  of  this  rcscar-h  effort  was  to  evaluate  the  six 
given  forecast  in, i;  methods  using  a-.,  actual  sample  i'OO  iloms  ttLookod 
at  a  base  level  consolidated  supp' .  a'livity.  Tlu-  selection  of  these 
particular  six  forecasting  nethods  was  based  on  choosin;.’,  methods  which 
could  be  used  routinely  and  'nexpensively  to  forecast  demand  for  a 
large  number  of  relatively  inexpensive  Items  based  entirely  on  hi.storical 
demand  data.  The  exponential,  smoothing  method  with  an  adaptive  smoothing 
constant  was  the  only  one  of  the  six  methods  evaluated  which  explicitly 
accounted  for  the  effects  of  exogenous  variable  (var iablcsother  than 
historical  demand),  yet  the  data  requirements  are  the  same  for  this  method 
as  for  the  other  five  methods.  Other  more  sophisticated  methods  which 
employ  such  techn iques  as  spectral  analysis,  iterative  dynamic  programming, 
Bayesian  analysis,  etc.,  were  rejected  as  being  impractical  for  base 
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levol  use  boc.iusi'  of  Llie  iiiy’her  Uvcls  of  conpulcr  stor.i;;i.-  .Tnd 
tional  rcqu irenenls . 

MovIur  Avcr.'iRos  Mctliod.  The  process  of  con  puling  .a  moving  aver¬ 
age  is  quite  straigiu  forward .  The  nalhetnaticnl  formula  for  calculating 
the  expected  demand  in  a  period  sinply  averages  the  demand  experienced 
in  tlio  previous  n  periods.  The  average  demand  iorecasti  minimizes  the 
sum  of  the  squares  of  the  differences  between  the  most  recent  n  obser¬ 
vations  and  i-he  expected  value  of  irhe  demand  in  the  period  being 
11  ' 

forecasted.  It  is  simple  and  easily  adapted  to  automatic  data 
processing  equipment.  Tlte  formula  is  as  follows.  ■ 


V 


t 


i:' 

i=t-n 


d. 

1 


;  (1) 


where:  v^.  is  the  forecast  of  demaml  in  period  t.  : 

I. 

d.  is  the  actual  demand  in  period  i. 

1  > 

n  Is  the  number  of  periods  of  actu.al  demand  which  arc  used 
to  develop  the  forecast. 


For  example,  let  the  period  be  one  month  and  n  =  10.  Then 

for  the  first  month  after  time  t  =  0  is  v.,  and  —  the  sun 

1  10 

in  months  -9  through  0,  i.e., 

0 


V 


1 


10  t=y  i 


i(d_g  +  d_g  + 


•  V 


the  forecast 
of  the  demands 


(2) 
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Tlte  effect  of  moving  average  forecasting  dt  pemls  on  the  I'.iir.ber  of 
periods  (n)  used  to  make  tl»e  forecast.  When  a  large  number  of  periods 
are  used,  the  weight  given  to  eacl»  period  is  relatively  small,  and  ran¬ 
dom  Fluctuations  have  little  effect  on  the  forecasts.  The  moving 
.average  metViod  Ls  more  sensitive  to  changes  in  demand  if  n  is  small. 

Tile  moving  average  metliod  works  best  when  the  demand  generating  process 
is  stable.  If  tlie  demand  generating  process  has  a  constant  mean  and 
variance,  and  successive  observations  are  uncorrelated,  tlien  the  sample 
moan  is  an  unbiased  estimate  of  the  process  mean.  If  these  assumptions 
are  net,  then  the  demand  generating  process  can  be  considered  stable. 

Least  Squares  Method.  The  least  squares  method  of  forecasting 
assumes  that  the  demand  dj,  in  period  t  can  be  predicted  using  w^.  = 
a  +  bt,  where  a,  b  arc  determined  from  historical  data  by  minimizing 

F  =  Z  (d  -  w  )2  =  Z  (d  -  a  -bt)2 

t  t  t  c  t 

In  the  above  equation,  w^  is  the  forecast  demand  Cor  period  t  and  d^. 
is  tlie  actual  demand  for  period  t.  Assume  that  a,  b  are  to  be  deter¬ 
mined  by  using  the  demands  in  the  previous  n  periods.  The  current 
time  t*  is  a  review  time,  and  that  the  time  period  from  t*  to  t*  -  T, 

T  being  the  time  between  reviews,  will  be  referred  to  as  period  0,  the 

period  from  t*  -  T  to  t*-  2T  as  period  -  1,  etc.  Thus  the  number  of 
the  n  periods  to  be  iir.ed  in  determining  a  and  b  will  be  -(n-1),  ..... 
-1,  0.  To  determine  a  and  b,  set  3F/3a  =  3F/3b  =  0  and  solve  the 

♦ 

resulting  pqnationr.  for  a  and  b. 


-  (n-1 ) 

9F  =  -2  ^  i(d  -  a  -  bt)  =  0 

7a  t=0  ^ 


(4) 


-(n-1) 

3F  (d  -  a  -  bt)  =  0 


db 


t=0 


(5) 


(n-1) 


-(n-1) 


where:  Z  »  -n  ( n- 1 ) ;  ^ 

t=0  2 


]_  n(n-l)  (2n-l) 
6 


and 


-(n-1)  -(n-1) 

Ua^  d  ;P=X 


c=0 


t=0 


(6) 


Then  the  solution  to  the ’ equations  (4)  and  (5)  are: 


12 


so  that 


n(n-l) (n+1) 


w  =  U  +  h  1 1  +  n-T| 

n  L  2j 


.  [P  +  rM.  U|  ;  a  =  U  +  b(n-l) 

.•  L  2  J  n  2 


w  =  U  +  12 
t  r 


p  +  n-1. 
2 


n(n-l) (n+1) 


t  +  n^  ,  t  =  1,2,3,.. 


(7) 


Equation  (7)  allows  computation  of  the  forecasted  demand  for  period  t, 
t  =  1,2,3,...,  using  data  for  the  demands  in  periods  -(n-1),  ...,  -1,0. 
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For  example,  let  the  period  be  one  month  and  n=10.  Then  the 
forecast  for  the  first  month  after  time  t*  =  0  is  and  is: 


w  -  U_  +  12 
^  10 


!  iu 


+  2  It  +  ii 
990  JL  Zl 


-(n-1)  :  -(n-1) 

where:  U  d  ;  P  “2  td 


t-0 


t=0 
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The  least  squares  metliod,  sonetinie.s  called  simple  linear  regression 


incorporates  the  ability  to  follow  a  demand  pat  ;ern  which  has  a  linear 
trend.  The  demand  process,  as  taken  fron  the  above  derivation  is  assumed 
to  be : 

=  a  +  be 

The  demand  in  period  t  Is  a  function  of  a,  the  intercept  value  of  the 

trend  line;  b,  the  slope  of  the  trend  line  and  t  the  number  of  periods 

"k  14 

between  the  ordinate  and  time  t  . 

It  is  possible  to  have  least  squares  metliod  forecasts  that  are  less 
than  zero  using  the  given  formulas.  It  was  tiierefare  reasonable  to  impose 
a  requirement  that  demand  forecast  must  always  be  greater  than  or  equal  to 
zero.  If  the  forecast  demand  (using  the  formula)  was  less  than  zero  for  a 
particular  period,  ttie  forecast  demand  was  set  to  zero.  This  procedure  also 
had  important  consequences  for  subsequent  forecast  demands,  since  the  fore¬ 
cast  demand  that  v/as  set  to  zero  was  used  for  forecasting  demands  in  the 
following  n  periods. 

The  moving  average  and  least  squares  methods  for  forecasting  demand 
both  require  demand  data  be  available  for  n  periods  back.  This  dis¬ 
advantage  is  not  inherent  in  the  exponential  smoothing  methods  to  follow. 

Single  Exponential  Smoothing  Method.  Single  exponential  smoothing 
is  a  forecasting  method  similar  to  moving  averages  in  the  sense  that  it 
provides  accurate  forecasts  for  a  stable  demand  pattern.  The  forecasting 
formula  Is: 

^t  “  "''t-l  ■  “^^t-1  ’ 

where:  =»  the  single  exponential  smoothing  forecast  for  demand  in 

period  t, 
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dj.  j  =  l  lie  .-lotiiil  clomiiul  i  X()or  icmicd  ;n  period  t-1, 

2^  =  tile  single  exponeiii  ial  smoothirg  forecast  for  deioand  oace 
for  period  t-1. 

a  =  the  sinoorhing  constant,  0  ^  J 

Tile  effect'iof  singl>*  exponential  smoothing  depends  on  the  size  of 

a,  the  smoothing,  constant.  If  a  is  large,  nore  weight  is  given  to  the 

most  recent  demand  experience  and  the  forecast  is  more,  sensitive  to 

fluctuations.  Instead  of  weighting  each  past  observation  equally  as  does 

■< 

the  moving  average  method,  the  weight  given  to  previous  observations 
decreases  geometrically  with  age.  If  the  demand  generating  process  has  a 
constant  mean  and  variance,  and  successive  observations  arc  independent, 
then  tlie  expectation  of  the  single  exponential  forecast  is  equal  to  the 
expectation  of  tlic  mean  of  the  demand  generating  process.  If  the  assumptions 
about  tlio  demand  generating  process  are  met,  single  exponential  smootliing  can 
be  used  to  produce  an  accurate  estimate  of  demand.  Therefore,  single  exponen 
tial  smoothing  is  as  accurate  as  the  moving  averages  technique,  while  coaputa 
tion  is  simpler  and  storage  requirements  are  reduced. 

Uouhle  Exponential  Smoothing  Method.  Oouble  exponential  smoothing  is 
a  forecasting  method  which  uses  a  first  order  polynomial.  Thus  it  incorpo¬ 
rates  the  ability  to  follow  a  demand  pattern  v/hich  has  a  linear  trend.  The 
demand  generating  process  is  assumed  to  be; 
dj.  =  a  h  bt. 

The  demand  in  period  t  is  a  function  of  a,  b  and  t.  TTic  coefficient  a  is 
the  Intercept  value  of  the  trend  line;  b,  the  first  derivative  of  d^  with 
respect  to  t ,  is  the  slope  of  the  trend  line;  and  t  is  the  number  of  periods 
between  the  ordinate  and  period  t*.^^ 


In  order  Lo  devc'Jnp  ilu'  IDree.isL 
tlal  smoothinv;,  coefficients  a  and  b  can 
statist  ics. 

®t  "  ‘‘t-1 


b  =■ 

t 


b 

c-1 


+  a^(Y  -d  ) 
t-l  t-1 


for  tb 
bo  est 


s  process  usLii^  exponca- 
.mated  using  two  smoothed 


The  one  period  double  exponential  sniootliing  for.jcast  is  then; 


Y  -  a  +  b 
t  t  t 


(9) 


The  forecast  for  t  periods  in  the  future  is: 


+  bj.T.  (10) 

VJrictcn  in  terms  of  single  exponential  smoothing  values  and  a,  the 
double  exponential  smoothing  forecast  for  period  t  is: 

Yj.  =  a  Xt  +  (1  -  (H) 

Use  of  the  double  exponential  smoothing  model  should  be  based  on  cer¬ 
tain  assumptions  about  tlie  demand  pattern.  Tf  the  time  scries  has  a  trending 
average  demand  rate,  either  increasing  or  decreasing,  the  double  exponential 
smoothin;,  method  produces  an  accurate  estimate  of  demand.  For  such  a  demand 
pattern,  the  moving  averages  and  single  exponential  smoothing  methods  lag 
behind  a  trend,  while  the  double  exponential  smoothing  method  is  more 
responsive  in  that  it  accounts  for  a  trend  factor, 

Triple  r.xponcnt  ial  Smoothing  Method.  Triple  exponential  smoothing 
is  based  upon  a  second  order  polynomial  time  series  process  which  ac¬ 
counts  for  tron<l  and  the  rate  of  change  of  trend.  The  assumed  demand 
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generating  process  is: 


=  a  +  bt  =  ^}  cf^  ; 

where  b  is  the  first  derivative  of  with  respect  to  t  evaluated  at 
t»0;  and  c  is  the  second  derivative  of  d^  with  respect  to  t,  which 
makes  it  the  coefficient  for  the  rate  of  change  of  trend.  Coefficients 
a,  b  and  c  can  be  estimated  by  computing  tliroe  smoothed  statistics: 

^t  = ‘^t-1  ' 


b(.  =  +  c  ,  -  ^^(2  -ot)(Z  -d  ) 

^  ^  t-1  2  t-1  t-1 


c^  =  c 


t-i  - 


where  ^  is  the  triple  exponential  smoothing  forecast  made  for  period 
t-l.  The  forecast  for  period  t=l  is  then: 


Z  =a^+b  +^0. 

t  t  t  t 


The  forecast  made  for  r  periods  in  the  future  is: 


Z^  =  +  b^T  +  h  c^t2  (13) 

Written  in  terms  of  double  exponential  smoothing  values  and  «,  the 
triple  exponential  smoothing  forecast  for  period  t  is: 

^t  “  +  (1  -  Vi 

Triple  exponential  smoothing  is  accurate  when  the  process  to  be 
forecasted  can  be  adequately  represented  as  a  quadratic  function  of 
time.  Thus,  if  there  is  reasonable  justification  of  assuming  that 
the  time  scries  can  be  represented  by  a  polynomial  of  the  form 
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r 


y 


=  .1  (-  M  +  ^  ct^, 

Clion  triple  cxpononti.il  .smootli inn  provide.*? 
20 

forecast  methods. 


tl^e  iio.st  .iccuratc  and  stable 


NOTE :  Eor  the  niven  three  .ind  higher  orde  •  oxnoiientlal 
snoothlnn  niotliofls,  Dro?;n,  in  his  proof  of  lie  runda.Tcnt.il 
theorem  of  exponenti.il  smootiiino,,  ha.s  denoistrated  that 
It  is  possible  to  estimate  the  (nl-1)  coef f  icients 'In  an  nth 
order  polynomial  for  the  corresponding  order  of  exponen- 
ti.nl  smoothing;  hy  usinp,  linear  combinations  of  the  first  (n-:-l) 
orders  of  exponential  snoothins.  It  shonlti  be  noted  .nt 
thi.s  point  tli.nt  brown's  proof  c.nn  be  used  to  derive 
formula.s  (8),  (11),  and  (14).^^ 


Adaptive  Ex n^o nential  1 1 inp. .  ;  T1 1 c  s i n p  1  e  exponential 

smoothlnp  forecasting  cpu.ntion  is: 

\  =  ^‘It-l  +  Cl  -  Xfl  C15) 


Let  X*  equal  the  adaptive  single  exponcnti.il;  smoothing;  forecast  for  period 
t,  and  let  a*  be  an  adaptive  .smoothing  consthnt.  The  adaptive  single 
exponential  smoothing  formula  is  therefore: 

I 

I  * 

X*  =  +  (1  -  «■*)  ^t-1  C16) 

I  ’ 

The  single  exponential  smoothing  method  provides  good  forecasts 
wlicn  the  demand  process  has  a  constant  mean  and  variance.  However, 


k 


J 


actual  demand  is  often  :iul)iect  to  sudden  inttemental  changes  due  to 
variables  exo>>enous  to  the  model.:  For  example,  the  number  of  systems 
for  whicli  an  item  is  used  may  be  suddenly  increased  or  decreased.  Under 
these  conditions,  as  the  level  of  demand  changes,  the  single  exponential 
smoothing  formula  will  change  the  forecasts  over  tine,  moving  towards  the 
new  demand  curve  at  a  rate  dependent  on  the  value  of  a,  the  smoothing 
constant.  The  adaptive,  aspcct.s  of  tlii.s  adaptive  single  exponential 
smoothinr,  method  were  built  upon  tliis  feature. 

The  basic  philosophy  was  to  change  a*  to  a  value  of  one  when  sudden 
shifts  in  the  actual  demand  level,  d,  were  detected.  Thi.s  has  the 
advantage  of  changing  the  estimates  of  the  mean  demand,  X*,  almo.st 
immediately  to  the  range  of  the  new  demand  mean.  The  procedure  to 
accomplish  this  was  based  on  tlie  observation  of  actual'  demand  values 
outside  an  acceptable  range  around  the  current  forecast. 

It  is  assumed  for  the  purposes  of  this  model  that  the  random 
fluctuations  about  the  basic  demand  level,  d,  were  normally  distributed 
with  a  known,  constant  standard  deviation.  In  such  a  case,  it  is 
natural  to  use  a  set  of  controls  on  X*  which  permit  a  rapid  response  to 
abrupt  changes  in  demand  patterns.  There  is,  however,  a  difference 
between  the  problem  here  and  the  apparent  analogy  with  the  classical 
quality  control  case  of'  controlling  a  process  by  detecting  points  out¬ 
side  control  limits  sec  around  the  desired  process  mean.  The  difference 
is  brought  about  because  the  initial  demand  forecast  is  likely  to  be 


difft'roiit  from  tlie  Lrue  moan  domaiid  by  an  amount  ilepondont  on  tao  'jc.ooth- 
ing  constant  and  tlio  demand  seqiicnce.  This  possible  error  means  that  a 
single  actual  demand,  d^,  may  differ  sufficiently  front  the  current  fore¬ 
cast  to  lead  one  to  suspect  that  tlie  basic  demand  has  shiiced  when»  in 
fact.  It  has  not. 

Because  of  this  complication,  two  criteria  were  developed  for 
concluding  that  a  new  demand  level  had  been,  detected  .vid  clianging  the 
smoothing  constant  is  necessary.  The  first’ of  tliese  criteria  is  ap¬ 
plied  in  the  case  v.Iierc  a  single  demand  is  discovered  outside  a 
range  of  four  standard  deviations  around  the  current  r\e.an  of  actual 
demand.  The  second  criteria  is  applied  when  two  successive  demands 

I 

differ  from  the  current  estimate  iiy  more  than  1.2  standard  deviations 
and  both  of  these  points  are  eitticr  above  or  below  th^  mean  actual 
demand. 

The  second  consideration  in  tiie  development  of  the  adaptive 
aspects  of  the  model  was  the  choice  of  the  values  of  tlic  smootliing 
constant  after  the  detection  of  a  shift  in  the  basic  demand  level. 

The  model  uses  a  3-stage  cliange  of  the  smoothing  constant .  The 
constant  is  set  equal  to  1  for  the  first  period  after:  a  change  In  the 
basic  level  of  demand  has  been  detected,  and  to  .8  Cor  the  second  period 
and  .5  for  the  third  period.  In  tlie  fourth  period  the  value  of  a  is 
returned  to  its  "normal  level.  ‘  I 


Clioosing  the  Smoothing  Constant 
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Approach.  The  amoothinR  constant,  A,  is  used  in  eac!i  of  cue 
exponential  snioothlnR  forecasting  methods.  According  to  Brown,  the 
value  of  ot  is  usually  chosen  bctvieen  .01  and  0.3.  Values  between 
0.3  and  1  are  normally  not  used,  since  these  values  cause  the  formula 
to^be  too  sensictye  to  random,  non-representative  fluctuations  in  the 
demand  generating  process.  For  the  exponential  tmoothiug  methods 
evaluated  in  this  study,  values  of  ot  between  .05  and  .5,  in  increments 
of  .05,  were  considered.^  The  value  of  et  which  minimized  the  standard 
deviation  of  forecasting  errors  for  the  single  exponential  method  was 
chosen  as  the  optimum  A. 

The  purpose  of  the  smoothing  constant  is  consistent?  in  all  of 
the  exponential  smoothing  methods.  This  purpose  is  to  weight  the  past 
demands  used  in  the  forecasting  formula  In  a  single,  prescribed 
manner.  Tlie  chosen  should  strike  a  balance  between  responsiveness 
to  demand  process  changes  (value  of  ei  approaching  0.5)  and  damping  of 
random  fluctuations  in  the  demand  process  (value  of  Ot  approaching  .05). 
It  sliould  also  be  noted  that  in  this  study,  the  same  set  of  demand  data 
was  operated  on  by  all  forecasting  metliods.  Considering  the  purpose 
of  ot  and  the  use  of  the  same  demand  data,  it  was  reasonably  assumed 
that  the  optimum  value  of  ot  for  single  exponential  smoothing  approxi- 
oiates  the  optimum  value  of  ct  for  double  and  triple  exponential  smoothing 
The  ot  that  minimized  the  standard  deviation  of  forecast; errors  for  the 
single  exponential  smoothing  method  was  therefore  also  used  as  the 
optimum  ot  for  the  double  and  triple  exponential  smoothing  methods. 
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This  optimum  «•  wns  also  usod  as  the  "normal"  value  of  A  for  tne  adap¬ 
tive  triple  exponential  smoothing  method. 

Method.  Specifically,  the  optimum  A  from  the  set  '.05  £  OH  0.5 
in  increments  of  .05,  was  selected  in  the  following  manner: 


(1)  X^‘  =  +  (1  -  oO  X^_^, 


(17) 


where  X  =  the  single  exponential  smoothing  forecast  for  period  t. 

t  i 

Xj.  =•  single  exponential  smoothing  forecast  for  period‘t-1,  except 

K  where  X„  =  —  -51  d  . 

t-1  0’  0  10'  i— 9  i 


(2)  Find  for  each  item  for  each,  value  of  A,  .05  <  CH  ±  0.5. 
There  are  10  values  of  X^,  one  for  each  o(,  for  each  itei®  per  month. 

This  was  repeated  for  all  316  items  for  each  month  of  the  forecast 
period,  1  November  1973  through  31  March  1974. 

(3)  The  differences  between  actual  demands  d^  and  forecast 
demands  X^  for  each  month  of  the  forecast  period  were  found  for  each 
item,  for  each  0^- 


Let  0(.  for  j  =  1,2,3, .,.,10  and 

®fly  °  \  ■  ‘^i*’  ”  1,2,3,. ,.,10) 


(13) 


Repeat  for  each  item. 


(4)  The  mean  for  each  e  was  then  found. 

J  i 

Mi  -z  .  n  =  312  items.  Repeat  for  each  (X  . 

J  n  i=l  WJ  1  j 


(19) 


(5)  The  standard  deviation  for  each  e„  was  then  calculated. 
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the  initial  coefficient  was  fictermined  by  taking  an  average  of  the 


actual  demand  over  the  10  month  base  period. 

0 


^0  ’  ll  \ 
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'  Initial  Estimate:  Doublo  Exponential  Smoothing,  the  double 
exponential  smoothing  method  requires  initial  estimates  tor  both  the 

ii 

slope  and  the  intercept  coefficients  a  and  b; 

Yj.  =  a  +  bt  ■  (2; 

Hie  method  used  to  compute  was  a  simple  linear  least  squares 
regression  based  on  the  10  data  points  of  the  base  period.  The 
least  squares  estimates  of  a  and  b  served  as  the  components  of 
for  the  forecast  The  remainder  of  the  5  months  of  Che  fore¬ 

casting  period  were  forecast  using  the  standard  double  exponential 
smoothing  formula. 

Initial  Estimate:  Triple  Exponential  Smoothingi  The  triple 
exponential  smoothing  method  requires  initial  estimates  for  the 
intercept,  slope  and  rate  of  change  of  slop  of  the  demand  process 

2K 

equation.*-'  Consider  the  assumed  demand  process: 

*  a  +  bt  +  *5  ct^  (2; 

Over  the  10  data  points  of  the  base  period, 

d  (demand)  ■ 

bt  +  ct^ 


-9  -8  -7  -6  -5  -4  -3  -2  -1  0  ' 


Yo 


Oemand  Versus  Time  Distribution 


The  value  of  a  is  the  d  axis  intercept  of  d  ,  and  can  be.  estimated 

t 


by: 


10  '’V 


(24) 


The^ value  of  b  is  the  limit  as  At  approaches  zero  of  — .  Using  data 
■  At 

points  at  t  =*  -9, -7, -6, -4,  =  3  and  -1; 


b  =  i 

3 


^  {U-U)  ^ 

- 2 -  - 2 - '■  - - 1 


J 


which  reduces  to: 


b  =  1  (d  -d  +d  -d  +d  -d  )  (25) 

6  -9  -7  “6  -4  _3  -1 


The  value  of  c  is  equal  to  the  rate  of  change  of  b  with  respect 


to  t.  Consider  Che  rate  of  change  of  b  using  two  data  points  of  b; 


c  .  1  1  i  |W-, 

1 

+  - 

7d_^-d_p  -(d.,-d_p- 

2 1  ^  r~^  2 

3 

2  2 

_  ^ 

4 

which  reduces  to: 

c  =  -i  (d  -d  -d  +d  )  (26) 

12  _9  -7  -3  -1 

The  following  equation  was  taken  from  Brown’s  book.  Smoothing. 

Forecasting  and  Prediction  of  Discrete  Time  Series,  and  is  the  triple 

exponential  smoothing  demand  process  quadratic  formula  in  terms  of  a, 

26 

b,  c  and  0(.  In  this  form,  the  demand  process  equation  can  be  used 

to  estimate  Z^,  F«r  Z  then; 
t  0 
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! 

t 


Z  =  a  -  3(1  -  «»)  b  +  3(1  -  «)(4  -  3a)  ,  (27) 

0  A  20?  . 

V  * 

I 

The  range  of  values  taken  on  by  must  be  1  mited  to  positive 
demands.  It  Is  also  reasonable  to  impose  an  uppe.-  limit  of  2a,  so 
that:  > 

0  _<  Zq  _<  2a 

When  Z^  calculated  by  formulas  (24)^  (25),  (26)  aad  (27)  is  negative, 
let  Zq  =  0.  When  Z^  calculated  by  ‘(24),  {25),'  (2f>)  and  (27)  is  2.  2a, 
let  2  =  2a. 

0  i 

The  value  of  Z^  at  t-1  =  0  was  calculated  in  the  manner  pre¬ 
scribed  above.  The  values  of  Z^  for  subsequent  periods  of  course, 
were  calculated  using  the  standard  triple  exponential  smoothing  formula. 

Evaluation  Criteria  For  Forecasting  Methods 

A  comparative  evaluation  of  the  accuracy  of  the  six  forecasting 
techniques  was  performed  in  terras  of  the  forecast  error  each  produced. 
Three  statistical  measures  of  error  were  selected:  the  mean  square 
error,  mean  forecast  error  and  variance  of  error. 

The  Mean  Square  Error.  The  mean  square  error  (MSE)  Is  a  measure 
of  forecast  accuracy  which  emphasizes  large  errors. 22  xhe  mean  square 
error  for  each  item  and  forecasting  technique  was  computed  by  summing 
the  squared  forecast  error  for  each'of  the  5  months  and  dividing  by  the 
number  of  months  forecast.  The  equation  for  the  MSZ  is  developed  from 
the  definition  of  forecast 

5 


error : 


t 
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where  o  =  the  absolute  avorape  forecast  error  ever  the  3  noath  forccastiii^ 

period  •  *  .  ; 

R  ■  V  ,  W  ,  X  ,  Y  ,  Z  ,  or  X*,  depending  on  the  forecast  sethexi 
t  c  t  t  t  t  • 

used.  I 

d^  =  the  ajctual  demand  in  period  t.  ^ 

Then  the  mean  square  error  is: 

2 


MSE 


(R  -  d J‘ 

t  t 


The  mean  square  error  for  each  item  was  summed  and  divided  by  tl:c  sample 
size  to  obtain  an  average  mean  square  error.  The  forecasting  technique 
with  the  smallest  average  mean  square  error  was  termed  the  most  accurate 
in  tliat  it  minimized  the  number  of  large  errors.^’ 

The  Mean  Forecast  Error.  This  statistic  was  computed  by  averaging 

t 

the  absolute  value  of  the  forecast  errors  fqr  each  item  by  forecast  model. 
The  equation  for  the  mean  forecast  error  is:’ 


\...  ' 


HFE  =  -L  e.  =*  -i-  jf  jR  -  d^l.  " 

n  t^l  j  n  t’n  '  t  t'  ; 


The  average  mean  forecast  error  i.jas  then  computed  for  eacli  forecasting 
technique  by  summing  the  riFC  for  all  items  and  dividing  by  the  number  of' 
Items.  The  preferred  model  was  that  for  which  the  average  Ml'E  was  closest 
to  zero. 

The  Variance.  The  third  comparative  statistic  was  the  variance  of 
forecast  error  (VAR).  A  high  variance  indicates  an  unstable  or  erratic 


estimator,  while  a  sm.ill  variance  Indicates  a  consisrcnc  est:na:or 


29 
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The  variance  of  forccnsit  error  was  computed  for  each  iteni  for  each' 
forecast  model.  The  equation  for  the  variancaaiy^  forecast  error  is: 

VAR  =  — -  y;"  (e  -  MFE)^ 

n-i  t 

s 

The  variances  for  the  sample  were  then  averaged  for  each  nodel.  Tlie 
model  with  the  smallest  average  variance  was  considered  the  most  stable 
est  imator . 

Statistical  Inference.  Ac  this  point,  it  is  pertinent  to  discuss 
statistical  inferences  concerning- tlie  error  statistics  of  the  six 
evaluated  forecasting  metliods.  A  nonparanetric  test  was  necessary  to 
evaluate  the  error  statistics  since  no  specific  assumptions  could  bo 
made  about  the  distribution  of  the  error  population.  Vonparametric  tests 
are  not  concerned  with  parameters  of  a  specific  type  of  distribution,  and 
no  assumptions  are  made  about  the  populations  sampled. 

The  nonparamotr ic  Kruskal-Wall is  test  (H  test)  v;as  used  to  determine 
if  there  was  a  significant  difference  in  the  effectiveness  of  the  six 
forecasting  methods  evaluated  in  this  study.  The  11  test  evaluates  the 
null  hypothesis  that  k  independent  random  sample  distributions  come  from 

identical  popula' ions  against  the  alternative  that  the  means  of  the 

30 

populations  are  not  all  equal. 

In  this  case,  the  absolute  average  error  distributions  generated  by 
each  forecasting  method  were  considered  the  independent  sample  distributions. 
Within  each  subset,  the  error  distributions  were  compared  to  determine  if 
they  did  or  did  not  come  from  Identical  populations.  If  the  pertinent 
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error  d  istr ibut ions  were  [  iwn  identical  popnlat Lons ,  then  there  were 
no  significant  differences  in  the  effectiveness  of  the  respective  fore¬ 
casting  methods.  On  t’.ic  other  hand,  if  they  were  not  from  identical 
polulations,  then  there  were  significant  dif f erc'nces  in  effectiveness. 

The  combined  hb.solute  average  errors  from  the  forecasting  methods 


being  compared  were  ranked  according  to  size.  was  the  sum  of  the  ranks 


assigned  to  the  observations  in  the  ich  method  and  n  =  + 

+  n  .  The.  test  v;as  based  on  the  statistic 

2 

1!  =  i  —  -  3  (n-kl) 

n(n+l)  nj^ 


If  the  null  hypothesis  was  true,  the  sampling  distribution  of  this  statistic 

could  bo  approximated  closely  with  a  chi-square  distribution  with  k-1 

degroe.s  of  freedom.  Consequently,  the  null  hypothesis  at  the  level  of 

2 

significance  7.  could  be  rejected  if  H  exceeds  "X  ^  for  k-1  degrees  of  freedom 
Aii')(^05  was  used  to  give  a  confidence  interval  of  95  percent. 

Final  Sample  Item  Cla.ssifications.  From  an  examination  of  the  toal 
demand  data  and  the  error  statistics  for  each  forecasting  method,  a  relation 
ship  was  noted  between  the  demand  cliaracteristics  of  an  ito.n  .^nd  the 
magnitude  of  its  absolute  average  forecast  error.  Items  with  high  eem.and 
rates  tended  to  liavc  larger  error  statistics  than  did  items  with  low 
demand  rates.  Consequently,  the  former  played  a  more  significant  role 
in  the  determination  of  the  best  forecast  model.  In  order  to  equalize 

^  .  V. 

the  effect  that  items  had  on  this  determination,  total  sample  of  items 
was  divided  into  three  subsets  based  on  Che  demand  proco.ss  cl'.araccerisc ics 
of  each  item.  The  six  demand  forecasting  methods  evaluated  were  therefore 
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applied  to  the  [uliowin.;  set  and  subsets  of  Che  316  icca  final  sac.ple. 

(1)  The  tot;il  set  of  31h  itens. 

(2)  The  subset  which  included  the  71  icons  with  relatively  ai^a 
level,  non-erratic  denand  processes.  Non-erratic  means,  in  this  study, 

a  relatively  constant  demand  process,  or  a  demand  with  a  discernible 

•jt 

t^end  (linear  or  curved)'. 

(3)  The  subset  w'hich  included  the  199  items  with  relatively  lov- 
levcl,  non-erratic  demand  processes. 

(4)  The  subset  which  included  the  46  items  wic'n  erratic  demand 
processes  (no  discernible  constant  level  or  trend,  linear  or  curved). 

Standardiy.at ion  of  Error  Statistics.  Consideration  was  given  to 
standardizing  the  absolute  average  error  statistics  used  as  a  basis  for 
evaluating  the  six  forecasting  methods.  For  example,  the  absolute  average 
error  for  each  Item  could  have  been  divided  by  the  average  monthly  demand 
rate  for  that  item  to  produce  a  standardized  statistic  independent  of 
demand  level.  However,  the  absolute  average  error  for  each  item  would 
have  been  divided  by  the  same  average  monthly  demand  rate  for  each  fore¬ 
casting  method.  The  standardization  process  would  not  change  the  relative 
ranking  of  the  forecasting  methods.  Standardization  of  the  absolute 
average  error  statistics  was  therefore  not  used. 
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(■rneral  Comparison  of  Foreca.;t 
Me t hod  Performance 


« 
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Frequenc 


0  .5  1.0 

1.5  2. 

0  2.5 

3.0  3.5,  4.0 

4.5 

Method 

Mean  Squared 
Error 

Mean  Forecast 
Error 

Variance 

- 

Single  Exponential 

Smoothing 

67.99 

5.51  ; 

38.15 

- 

Moving  Average 

76.18 

4.68 

44.49 

— 

Double  Exponential 

Smoo  th i ng 

103.19 

6.14 

66.38 

mm 

Least  Squares 

108.92 

6.24 

71.02 

- 

Triple  Exponential 

Smoothing 

124.67 

7.18 

74.21 

Adaptive  Single  Exponential 

120.78 

6.22 

83.28 

Frequency 


Subset  2  (199  Items) 

Absolute  Average  Error  Oistributlons- 


Average  Error 


0  .5  1.0 

1.5 

2.0  2.5 

3.0  3.5 

4.0  4. 

Method 

Mean  Squared 
Error 

Mean  Forecast 
Error 

Variance 

Single  Exponential 

Smoothing 

0.36 

0.49 

0.12 

Moving  Average 

0.37 

0.49 

0.12 

Double  Exponential 

Smoothing 

0.40 

0.50 

0.14 

Triple  Exponential 

Smoothing 

0.41 

0.49 

0.17 

Least  Squares 

0.47 

0.53 

0.20 

Adaptive  Single  Exponential 


0.52 


0.56 


0.20 


Double  Exponential  Smoothing 

3.56 

1.32 

1.86 

Moving  Average 

3.68 

1.31  ; 

2.01 

Single  Exponential  Smoothing 

3.81 

1.35 

2. 04 

Adaptive  Single  Exponential 

4.51 

1.51  • 

2.23 

Least  Squares 

4.46 

1.32 

2.7d 

:  5.52 


Triple  Exponential  Smoothing 


1.46 


3.47 


0  .5  1.0  1.5 

2.0  2.5 

3.0  3.5  4.0 

4.5 

Method 

Mean  Squares 
Error 

Mean  Forecast 
Error 

Vari ance 

- 

Single  Exponential  Smoothing 

16.02 

1.74 

13.02 

- 

Moving  Average 

17.84 

1.78 

14.72 

- 

Double  Exponential  Smoothing 

23.89 

1.89 

20.39 

- 

Least  Squares 

25.36 

1.92 

21.72 

— 

Adaptive  Single  Exoonential 

28.05 

1.97 

24,24 

Triple  Exponential  Smoothing 

29.00 

2.13 

24.55 
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IV.  Detailed  Summar 


Experiment  Summary 

Sampling  Tccbnlque.  There  are  many  possible  ways  of  seieccinj;  a 
sample  of  EOQ  ifems.  One  method  is  by  stock  niober  classification. 

Other  methods  are  demand  level  and  cost.  In  this  study,  the  strati¬ 
fied  base  level  X152  and  XB3  item  sample  (B  cost  group)  was  selected 
from  three  randomly  chosen  systems.  Within  the  toal  sample  set,  iteais 
were  grouped  according  to  demand  level  (subsets  1,  2  and  3).  These 
Important  facts  should  be  kept  in  mind  when  drawing  conclusions  from 
the  results  of  the  analysis. 

Evaluation  Criteria.  The  forecasting  methods  were  ranked  according 
to  error  statistics.  Variance,  mean  forecast  error  and  mean  squared 
error  were  the  error  statistics  used.  The  smaller  the  error  statistic, 
the  more  accurate  the  forecasting  method.  The  Kruskal-Wallis  test  was 
used  to  determine  whether  there  was  significant  differences  in  the 
relative  effectiveness  of  the  six  forecasting  methods. 

Findings 

Error  Statistic  Banking.  Single  exponential  smooching  proved  to  be 
the  most  accurate  alternative  to  the  present  method  of  moving  average 
forecasting.  Double  exponential  smoothing  also  ranked  in  the  cop  three 
for  the  total  sample  set  and  for  subsets  1,  2  and 

triple  exponential  smoothing  and  adaptive  single  exponent  iiicochi.'.g 
ranked  fourth,  fifth  and  sixth,  respectively.  The  following  paragraphs 
discuss  the  reasons  for  their  relatively  poor  showings. 


The  lon;>t  squares  netliod  tended  to  respond  to  sud  -ier. , 

non-representative  demands  by  excessive  changes  in  the  slope  of  the 

trend  line.  In  some  casfs,  this  cliange  of  slope  resulted  in  forecasts 

of  negative  demands  for  the  following  period  or  periods.  This  type 
£ 

of  erro^  was  mid'lgated  to  some  extent  by  setting  all  negative  least 
squares  forecasts  to  zero.  On  die  other  hand,  louble  exponential 
smoothing,  which  also  assumes  a  linear  trend  in  demand,  gave  relatively 
more  welglit  to  demand  history  with  subsequently  less  drastic  changes 
in  trend.  The  non-negative  demand  constraint  was  not  necessary  for 
double  exponential  smoothing. 

The  two  primary  causes  for  the  relatively  poor  showing  of  triple 
exponential  smoothing  were  the  absence  of  demand  processes  which  could 
be  accurately  approximated  by  a  second  order  polynomial,  and  the  poor 
results  obtained  with  the  formula  for  initial  conditions  None 

of  the  316  items  in  tlie  sample  had  a  demand  pattern  with  a  consistent 
rate  of  change  of  slope.  The  formula  used  to  calculate  initial  condi¬ 
tions,  though  mathematically  correct,  was  therefore  not  representative 
of  the  actual  demand  processes.  Starting  with  relatively  inaccurate 
initial  conditions,  the  errors  tended  to  compound  during  the  forecast 
period. 

The  results  obtained  with  the  adaptive  single  exponential  smooching 
were  also  disappointing.  The  major  cau.se  of  error  relative  to  standard 
single  exponential  smoothing  was  the  first  criteria  for  detactir.g  a 
change  in  demand  level;  i.e.,  .i  single  demand  detected  c.uts-Cc  four 
standard  deviations  of  the  mean  demand.  When  suc'i  a  ccaar.c  was  detec..eo, 
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the  smoothing  constant  viis  changed  to  one.  This  crtter^ia  and  the  subse¬ 
quent  smoothing  coustanc  ctiange  caused  the  Coretasc  to  change  excessively 
in  response  to  abrupt,  non-representative  demands.  Tlje-  second  criteria; 

i.e.,  two  successive  demands  differing  from  the  mean  demand  by  more 
* 

1,2  standard  deviations,  both  demantlsbeing  above  or  below  the  mean, 
appeared  to  change  the  smoothing  constant  at  appropriate  times  for  the 
given  sample.  .  i 

Three  obvious  approaciics  to  improving  the.  accuracy  of  adaptive 
single  exponential  smoothing  are  open.  The  first  would  be  to  drop  the 
flrst  criteria  for  changing  Che  smoothing  constant,  retaining  only  the 
second  criteria.  The  second  approach  would  be  to  change  the  raaxinun 
smoothing  constant  to  some  value  less  than  one.  The  third  approach 
would  be  a  combination  of  the  first  two  approaches.  However,  these 
possible  improvements  were  not  attempted  in  this  study  and  are  ieic  as 
subjects  for  further  research. 

Finally,  a  few  words  concerning  the  choice  of  smoothing  conitants 
for  double  and  triple  exponential  smoothing  are  In  order.  It  was 
assumed  that  the  optimum  smoothing  constant  for  single  exponential 
smoothing  would  also  be.  optimum  for  double  and  triple  exponential 
smooching.  No  proof  was  offered  for  this  assumption.  . It  would  be  of 
interest  to  compare  the  accuracy  of  double  and  triple  exponential 
smoothing  forecasts  using  smooching  constant  (o)  values  in  Increments 
of  .05  between  .05  <  a  ^  .5.  This  comparison  was  not  atteriptcb  for 
this  study,  and  is  also  left  as  a  subject  for  pertinent  future  researc.i. 
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Tlie  Kruskal -W.i  11 13  Tisi.  Tliis  u-st  indicaroil  that  the  populations 
of  absolute  average  foreirast  errors  are  nearly  [dentical  for  the  six 
forecast  methods.  liasod  on  this  test,  there  were  no  significant  dif¬ 
ferences  in  the  relative  effectiveness  of  the  s-fx  methods.  However,  as 
/• 

previously  no  d,  there  is  substantial  evidence  that  improvement  is 

required  in  tlie  formulation  of  le.ist  squares,  triple  exponential 

smoothing  and  adaptive  single  exponential  smootl.ing  forecasts.  It 

i 

would  seem  appropriate  ;tlu’n,  to  withhold  final  judgement  on  these  three 
forecasting  mechods  until  the  suggested  improvements  arc  made  and 
tested. 


Kocommendat ions 

The  purpose  of  this  resoar c!)  effort  was  to  examine  and  compare 
five  methods  of  base  level  demand  forecasting  which  ran  be  considered 
as  alternatives  to  the  present  metliod  of  unweighted  moving  averages. 

All  six  mechods,  including  unweighted  moving  averages,  were  evaluated 
according  to  relative  accuracy  in  describing  tlie  actual  base  level 
demand  patterns  of  a  statistically  chosen,  stratified  sample  of  XB2 
and  XB3  EOO  items. 

There  is  little  doubt  that  the  single  exporential  smoothing  method 
of  forecasting  item  demand  is  a  viable  alternative  to  the  present 
method  of  unweighted  moving  averages  for  items  represented  by  the 
sample  used  in  this  study.  Double  exponential  smooching  also  appears 

viable. 
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Research  uslnj;  the  given  forc.east  ing  methoiis  should  corcainly 
continue.  Rccoirmiend  that  subsequent  research  e.  forts  use  itea  s.iaples 
gathered  by  different  metho<lologies  at  base  level.  In  tiiis  way,  a 
broad  base  of  item  samples,  wholly  representat i’'e  of  all  EOO  items 

I 

Air  Force  wide,  can  be  established.  Given  that  the  results  of  future  \ 
research  efforts  bear  out  the  findings  of  this  study,  serious  consider¬ 
ation  ought  to  be  given  to  single  .and  double  ex;ioncnt  ial  smoothing 
methods  as  alternat ivos . to  the  present  method  of  forecasting  demands. 
Single  and  double  exponential  smoothing  have  the  well  documented, 
inherent  advantages  of  smaller  data  storage  requirements,  greater 
computational  ease,  and  greater  flexibility  in  adapting  to  demand 
pattern  trends  compared  to  unweighted  moving  averages. 

Further  research,  as  indicated  under  the  Findings  Section  of  this 
chapter,  is  necessary  to  properly  evaluate  all  pertinent  improvements 
concerning  least  squares,  triple  exponential  smoothing  and  adaptive 
single  exponential  smoothing.  Relatively  simple  adaptations  to  the 
FORTRAN  program  used  in  this  researcli  effort  can  be  made  to  analyze  the 
suggested  improvements.  Recommend  these  adaptions  be  made  and  the 
program  used  to  evaluate  the  given  sample  and  other  base  level  samples 
selected  by  different  metliodologies. 
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Appendix  A 

FO-tTRAN  Computer  Program 


-  A  - 


lOCNT  MAIN 

Flic  ?  »  luPitT  ^  IIMTT  a  - 

FILE  6  a  OUTPUT  P  UNIT  a  PRINTER 

_E4L€-LO-a_DSKF-IL»  -IJNlTL-a  llISK,WnRKElur.>RANtUjH,RECXIRX»_s_La^ - 

FILE  II  *  0SKrL2»UNIT=0ISK»WnRKrilE>RAN00H#RErCRDs3OO 

_ MMENSION  NSN<4>#QBSVO(15)»CALr.f  20)ApERICDUa)AS^0OTHC-SaXa - 

IC0NSNT(10),ST0RE(1^  ,HQUn<lO),ST\/(3)ASTW(3)»STX(3)ASTY<3)>ST2{3), 


B- 

C 

-C- 


alpha  nsn 

-frA4^^ZERO/a.O/^PLOAa/-l<3Aa/A4LPy<T.^/*KO«IHT/&/>NARKAtA 


INN  »  5 


ZEROINC-ABRAY-CAL-C  IN-  PREPARATION  FOR  CAUULAJEO  VALUED 


_ tm_2_tL=L»-20 _ 

CALC(N)  =  ZERO 
-2.JC0ttTlNUE. 


DO  3  M  s  IaAOO 

RTAPf  TMM.^^1  .TMn  a  A^  Ol  SN  (  T  3  «  t  a  I#  4.3  «  (  03  &3/I1LJ  3  A  J  *  Lt  1.5  ) 


51  F0RHAT(2X»AA4Al5r4,0) 

_ iCaUILt---  KOUfiT  -  1 - 


WRITE!  10>K0UNTa5?)  (MSN(  I),  I“  W4)A(0nSVn(  J)#  JsU  15)  A(CALC(K)AKaU 


52  F0RMATC4A4# l5F4.C»20r6.2) 

^  CQN'»’IMl!r 


4  HR[TE(inUT*60) 

-A4-F<}RKAT(IHU50X,31HSMn0THING  CONr^TAN-T  CAI.CUL ATI ONS.// A2-XA2UlNAT.tONA_ 
IL  STOCK  ‘NUMBER»7Xa5HA=.05»5Xa5HAs.10a5X>5HAs. 15a5Xa5HAs.20a5>  a5HA» 
- ?»25>5X/5HA  =  .30#^5X^5HAa»35»5X^5HA  =  »40#5X^5HAa»45»5X^5HA**50^AA) - 


00  5  Kal,lO 
-ALP-ALP^  r<>5- 


CONSNTCK)aALP 

-^-RIOOdOa-ZERO-- 


5  continue 

-24-BEAatl0^MARK.»S2)C14SNCU»I  =  l»4>*-C0asVDCJ.)*aaULSX»-CCALC£X)«Ka_LiL2iU — 

♦ 


Nsl 

^  A  M 


Xla(C0NSNT(K>*08SV0(Mn  +  ((l  .0*CnNSNT(K))*Vl) 

- ^OCTH(M)=Xl - 

DlFFs0irr+ABSCXl-Q9SV0(L)) 

- V-1-5.XI- - - 

NsN^i 

■cmwuc - - - -  ■  - 

CALC(K)=01FF/5.0 

^ER100(K>*PERIQ0<K)+CAI.C<K> - - - 

6  continue 

- 00  47-  IsU3 - 

STV(I)=0.0 

- STw(i)=a->o - - — - - 

STXCDsO.O  ; 

STZ(I)=0.0 

_ S.I2A(IJ  =  a.O  -  - , - 

47  continue 

61  FORMATdH  /3X»A^»3H  00  3  A  4  ^  6X/’r6 . 2^9  ( F  1 0 . 2  )  ) 

_ IW-tTEn0>MARKA52)CNSNn)W  =  1^4)^(00SV0CJ)^J=lA-,m^CtULUJtX.K=Uia3— 

NRlTE(ll2MARK#53)(SMnQTH( J)^ J=1^50) 

_54-WRHAT(50F6.2) - 

HARK=MARK+l 

- IE<.H.ARK.a.E.KnUN-I3  GO  TO  -21 - 

COUKTsKOUNT 

- 00-0-  I -1  MO  - - : - 

HOLCdJsZERO 

-0-SToRE(I)*PERIOO(r)/CnUNT .  -  - 

«RlTE(I0UT»62)(STQnE(M),Hst^l0) 


MARKai 

40««TlaC0UNT-l ,0  - 


aO  REAOClOiHARK»52)fNSN(!)»I=l»4)^(08SVO{ J),J=l,l5)#(CALC(K).K=l>2G) 

_ on  fll  .ls;t«LO - i _ 

OIFr»CALC( J)-STnRE( J) 

___EUUlE?=DirF-*DlFF.  . . . . . 

HOLOCJ)=HQLO(  J)*PLACF. 


MARK=MARK*1 

_ yrrManK.i  r^KniiMT^  an  Tn  no _ ^ _ 

HARKsl 

- WJ_e2-K  =  1 »  1 0- - - 

?LACE=SQRT(MQLD(K)/CnUNTn  • 

^ - H-CiQ-CiCJ-RLilCE - - - 

W  82  continue 

- MBlTECIflUI*4T)(Hni.QtH).^M=U10-) - 

63  FQRMAT(1HO»4X»20HSTANOARO  DEVIATION  10(4X/F6%2)>//> 

- PLACE =HOLO<-1) - ^ _ 

IP0INT=1 

■ - - - — 

IFtHOLDID.GE.Pl-ACE)  GO  TO  fl3 

- PLA&E-UOL-frT-I ) - - - 


IP0INT=1 

-  a^CONMNUE - - - 

ALPaCONSNK  IPOINT) 

- lia4TECinUT,64V- ALP. _ _ 

64  F0RKAT(1H0»44X,23H*-*«-**  ALPHA  SELECTED  s  ,F3.2>7H  ***••  ) 

_ KPOINT^O _ _ 

DO  84  K=l,lPniNT 

.84  KPXll-NI=KPaiHT»5 _ 

W  LP0INT=KP0INT-4 

_ liRJ.T£ClQUT/65)  .. _ 

65  FORHATC IHI,49X/ 35HSTNGLE  EXPONENTIAL  SMOOTHING  METHOD//) 

_ Nfl.IIEUnilT«66) _ _ 

66  F0RHAT(1H0/42X/14HF0RCAST  OEHAND/ 31X/ 1 3HaCTUAL  DEMAND/19X/ 

- LdHAasaLUTE) _ _ _ . 

WRlTE(lOUT/67) 

- 6JLXaRKAT-t  lH0/2X/2lHNATlaNAL-ST0CK  NUH8ER/ - 2tL0X^.7HN0V_iL__ 

17H0CC  /7HJAN  /7HFEfl  /7HHAR  )/2X/ 1 3HAVERAGE  ERRQR^//) 

- _ _ _ _ _ 

OIFFsO.O  ' 


,-3 


RCA0C10>MARX>52)(NSNa)»l3l>4)^(nBSVnCJ)/J=*l/l5>^£CALCCX)#K3l,2Q:, 

mm  ^  ^  m  M  IN  !<.«  .n>\\ 


JPQINTsLPOINT 

- 

L*N>10 

-4J4FF=0irFfA0S(SMC0THCJP&I»a)--0ft5.mCL-Xi - - - 

JPQINTsJPOINTf i 


OIFF=OirF/S.O 

_ mHTE(IQUI*&a)(NSN(I}>I=t^4UCSli0aiHXKXttC=Lg£iIH.'UXP0IUL: 

i£oasvo(j)/j=ii>i5)»oirF  » 

■■Aa-FilfiHA.T-ClH_j  2X,A^>  3H.  0Q_^ lA-AjL-diUSCr?. 2  1  QX*5  C  F7 ^ 
CALC(16)-0IFF 

CTvri  ^=«;TY/1 


STX£2)5STX(2)*01FF 

- 04FF-»  a*  0 - - 

VsO.O 

- M-IO-I-VaIO^ - 

ZEROING  ARRAY  CALC  IN  PREPARATION  FQR  CALCULATED  VALUES 


cqunt=i 

- P£R-IOD(  I  >*PLOAO  ♦COUNT- - • - 

to  continue 

- COUNT-KOUN-T - - - - - 

VlaV/lO.O 


00  11  J=2#5 

- 

L  =  J  +  9 

nasvatK-)-)  *nasvD.CL-: 

VlaV/lO.O 


ii  continue 


- M_t2  K=  I  #10 - 

Ua  U*Q05Vn(K) 

- |JLs^p*(PERI0D(K)*R8SVa(Kn - 

Pl*PU(PeRI00CR)*a8SVC(K*l ) ) 

_ pa*P2-*  4PE-izrQD<K)*afl4yctx*2U — 

P3aP3f(PERIQ0CK)*03SVn{K*3)) 

- P-AaP4-*-(PERICn(K)*08SV0CK*4V) - 

12  CONTINUE 


CAROS  NUMBEREO  UflS-  1725#  17/S5>  1805^  1845  PREVENT  LEAST 

— rRoH-omifl  neoative^  a-  negaiive  -i&.  set,  to-o^ - 

a?u/iQ.jQ _ ! _ 

Wla(CP*(CNNST*U)  3/990.0)aCl.04CfiNST)*l2*0  +  a  J 

■IFCHlAT^g^ai,  Hi=  (Uil _ 

CALC(6)=Wl 

iiau-oRsvnf  n*nBsmi.iti _ 

QaU/lO.O 

Jt2^CCPl-4CCNN5I«U>-U99a>a>aXl«.a*CtiNSTX*L2«4Ua - 

IFCH2.lt. 0«0)  W2=0.0 

jULC.Cn=H2 - - - 

(j8U-QBSVOC2)*aOSVOCI2) 


H33(CP2  +  (CNMST*tJ)  )/990.0)*C1.0*CNNST)tl2.0^Q  . 

.tfH.ViJ^U.a.-0V-H3=(Ua - 

CALC(a)=W3 

AT*y-O8SVO(3)*0nSVD(tVV - 1 

a=u/io.o 

^^4*((P3*(CNNST*lj->)y900>0)*^-l-.g^P-*<NST)*t;>>0»a — 
IFCM4.LT, 0.0)  H4  =  0.0 

CALC  ( 9  )=H4 _ - 

UaU-03SVOC4)*aOSVOC14)  ‘ 

^aU/XO.a _ 

H5aCCP4  +  CCNNST»U)  )/990.0)*(l. 0*CNNST)*1 2.0  +  0  '• 

■IF-Cli5.1  T.n.O)  H^gQ.n _ 

CALCC10)aW5 


j.  L«i  aU  a-C  QESlCaC  J-i  *.  J.5  L»i! 


85  COKTifiUC  - 

<iTyf9^=«-.T1tC?Wf.QU.tll - J - 

0IFF«0.0 

— inuT»6<j) - - - : - - — 

69  rORHAT(lHl,55X,21HHQ\/IHG  AVERAGE  METHOD^/) 

- t*fH-TE(lOUT*66) - 1 - 

HfliTEC 10UT»67)  “ 

- DO  e/  MARjUUKCUm^ - - 

REA0(  105MaRK,.S2)<NSN(  n»I  =  lA4)/(nBSVn(  J)»  J  =  l»15)»CCALC(K)>Ksl#2a) 

- S?€AO(  t  l>MAftK>5  3)-tSWOaTH(  i - 

00  83  N=l»5 


10- 


0IFr  =  DirF  +  A83(CAI.C(N)-n9SVn(Ln 
-COhT 


OIFFsOlFF/5.0 

-W41)sSTVn)^(niFF*niFFV - 

WRITE(  10UT,4a)(NSNCI),I  =  l,4)#(r.Al.C(K),K=l,5)# 

W^B5V0(  J)aJ=1  1>15)>0IFF - 

STV(2)sSTV(2)>0irF 

■  CALCn7VsO-ltl|: - 


INOsl 

- OiFF-aa^a. 


X«l. 0-ALP 
YQaCALCCai 


00 

-1 


89  LsLPQINT^KPOINT 
LP-nLSHflOTH.a  l.lt(X«Y01 


CALC(INO)aYl 
YQ^Yl - 


IN0=IN0+1 
-M-COUTINUr - 


WRITE(  10>MARK#52)<NSN(n^l  =  l»A)^(OBSV0(J)#Jsl>15)^(CALC(X)#K=l#20) 

-gy-CO^TPWE - - - 

STV(2)aSTV(2)/C0«iNT 

- - 


HRITECIOUTWO) 

-70-MflMAT(  lHU58XA2flHLEASl^  SOOARES-ME-tHOlW-)- 


HRITE(inUT>6fi) 

.;Aft4TE-L-l-nUT^-7-)- 


00  90  MARKsl^KoUflT 

-R€A0Cl0>MARK^52)(NSN<I)/^Is^l^4)»^<0WV0(J->^»lr#-l-5>»^XCAL4:<-X-^X^U:»O^ 


« 

1 

00  91  N=6/10 

- LaN*S - 

OIFF=DirF*ABS(CALC<N)-neSVOa)) 
gl.XOMTlHUC _ 


OIFFsDlFF/S.O 

lQUT^6aKNSNC.t)/t=!i^4-l/  CCAUCIiU  UlXt- 


1(QBSV0( J)/ J=1 15 )»DIFF 

-STwcz)-.st;u?)»otfe _ 


CALC(ia)=oiFr 

-SiW-C-U^  S.T  H  C  Li>  C  QI  Fij*  fl  I EXX 
DIFFaO.O 

- — - 

00  92  1=1/10 

-\taV*0!3S.\/0C-l) - - 


92 


continue 

-lOx-U/lO^O- 


AaVl 

-B«O&SV0  (  ll'-OOSVO^  XV+aB-SWU->-0W>i/0C6X*0flSAmCX4^0BSV4C-9->- 

6=8/6. 0 

Caoosvp(i>-QBsvo(a)«npsvn{7)-»QBSvPU) - 

CsC/l2.0 
-AtPHaALP - 


ALPH2=ALPH*ALPH 

-3WCi.*-X)/ALPK)«a _ 

Cls<(3.*X)*C4.-3.*ALPH)/(2.*ALPH2))*C 

-lOiAiflOXC-J _ 


IF<ZO.LE.O.O)ZO=C.O 
XEX2Q .. OT .. C2...0  * XI ). .  2 6=  2.>a*A- 


ZsALPH*CALCCn*(X*ZO) 

j[lALCCll.y=Z _ 

DO  93  J=2/5 

■Lajiin - 


ZTs(ALPIUCALC(  J))*(X*Z) 

jutc.aj=zi _ 


ZaZT 

-CON-TINUE  — 


MRrTE(lO>HAnK/52)<NSN(I)/I=t»A)/(nBSVo(J)/ Jsl/15)/(CALCCK)/Ksl/ij) 

90  continue _ _ _ 

STH(2)sSTH(2)/C0UNT 


J 


- I  ) - 

WR1TE(I01)T,67) 

- MARKaWKQUKT - 

REA0C10>MARK,52)(NSN(I)^Isl,4),(OflSVO(J), J«l,15)#<CALCCK)#Kal#20) 


LsN+lO 

_ 04-rr=airF>AQSccAUXfiX-4iflS.vaCL)-i- _ : - 

95  CONTINKC  • 

_ D^XF-ailirF/S.Q - - 

HRlTE(lOUT*6fl)(NSN(  n,I  =  l,4)#(FAI.C(K),Kal>5)#  , 


CALC(19)=DIFF 

_ SrUa-i-S-TYCD^COlFF^DIFF) _ . _ 

STY(2)=STY(2)+0IFr 

- auFsa.o - — - 

WRlTE(lO>MARK»52)(NSN(nM  =  l»4)#(OaSVO(./)» Jal»tS)^(CALC(K)»K«l>20) 


STYC2)sSTY(2)/C0l!NT 

- WR-imiOU-T-,-??).- - 

72  FQRMAT(1H1>4<>X#35HTRIPLE  EXPONENTIAL  SMOOTHING  METHOD) 

- V(«-ITE(-I0UT,66V - - 

WRITE(I0UT,67) 


REAC(lO>MARK/52)(NSN(n/Ial/A)^(08SVD(  J)»JsUl5)#(CALC(K>»Kal»2fl> 

-00-9-7  Nal  1^15  -  - - 

OIrraOIFF+ABS<CALC(N)-nBSVO(N)) 

XONTIUUC  - - - 

OIFFsOIFF/5.0 


U08SV0(  J)»J=1 1/15)»0IFF 

_ Cja.CC2C.)a0lFF_.. - ^ _ 

STZCDaSTZd  )KDlFF*niFF) 

_ ST2C2)aSTZ.{2)»Dirr _ 

HRlTEC10>MARK,52)<NSN<I)d  =  lr4)/(0BSV0(J)^J=ld5)#(CALC(K).Ksl>20l 
nTrr=n*n 


96  CONTINUE 

_ S-T2X2)sST2(2XyC0UN.T, _ _ _ 


WRlTE(inUT,73) 


_>4«H4EC40U^^46-V - - - - - - - - 

WRITECI0UT#74) 

^--P0«MAUlH0^^1HKm0NAL  -STOCK- NUf!BE4?^4X-^HA - WOV^SX^TMA— --OEC^^X^ 

17HA  JAN>5X^7HA  FEB^5X^7HA  HAR>5X>3WN0V# 4X> 3HDEC^ AX# 3H JAN* 4X# 


00  98  MARKsl/KQUKT 
— r^ccio^mark>-52 )<nsnci>».ixi #^A3>-CntLSunxj.> ^is | ,»is?#^crALrc"">>t<a^  r?^? 
REAOUl>MARK#S3)^SHOnTH(J)>Jsi>50)  ^ 

_4ua^a - 

00  99  1*1*10 

>..vsv»oflsvo(n _ _  - _ 

VlsV/lO.O  I  " 

--JCOaVl _ i 

AsALPH  ~~  ^  ^ - 


LslO 


00  100  J*K«L 

- D4fx=  oas  V  ac  4  x-- « t _ _ _ 

lOO  PLACE  =  PLACE-KOlFr*orrr) 

- ST0EVsSQRT(PLACE/9^)-. _ _ _ 

T£STsA,0*STDEV 

- M-2aVl.*TCST _ 1__ _ 

V3aVl-TCST 

- IfLCOaSVXUU*  t  X*GE«-V2  X  Aal  ^ _ _ _ 

IF<C8SVD(LM).LE.V3)  A  =  l.O 

- tF-tX»Ea.-V^O)  GO-  TO- 103 _ _  . _ _ 

IFCL.GE.IA)  GO  TO  10.3 

- TEST2-1  .g^STOC-V - ^ _ 

V2aVl+TCST2 

- T  E  S  T3- - - - 

IF(  (0BSVD(L+n,GE.V2).ANn.((]BSV0(L>2).CE.V2))  A*1.0 

- l£C.A*E0.-Uai  GO-TCL  tOl _ _ _ _ 

IF((0BSVDCL*1 ),LE.V3).AN0.(nBSV0(L+2),LE.V3))  4*1.0 

t03  Xal.O«A _ _ _ _ _ 

XTa(A*GnSVO<Ll)*(X»XO) “ 


-9 


STQREC2*K)»XT 
JtflaX-T - 


Va(V-OBSVO(Kn»ORSVO(L>l  1 

-Vlavyio.o _ 


ir(ST0RG(2*K-l).f:(}.1.0)  /.a.fl 
-lXCSIflREC2*K-l)-.Ea.*a)_A3»5 
lF(ST0RE(2*K-n.Ea..5J  A^ALPH 


IF  YOU  WANT  ALPHA  TO  SKIP  .8  WALUC*  CHANGE  THE  .fl  VALUE  CH 
XAfia-NUMaC  RE0~a5Sl^Tfl_A_  .^AIUUREM&iLE_ 


c 

w 


LaL*l 
■  IFa^LC-^ 


.M-)-GO-:tO-Sg- 


DIFFaO.O 

_H4..  t  L - 

DO  lOA  H=2j10»2 

F  aO I FF  ♦  AB  S  (  ST C R E  (  H4-OR S VO  C M-l-W — 


4^ 


QIFFsOlFF/5.0 

JtJt4JEClOUT*-7S)(NSN(n»laUAl«(S.TTlR£.CK.t»<a.L^lQl  .« _ 

ICOBSVOC J)^J=l 1>15)»01FF 

-TlS-TOftM  A-t(2X>A4>3H.00^SA4>10CF6*2  XjLjX^  SIF  T-^^X^-SXjX&^AX 
STZA(naSTZA(l)  +  (OlFF*niFF) 

- STZAl2  3sSUAf.2)*niFF - 


SH00TH(l)a01FF 

-lifi-ITECn2MARK^5S)-(SHnOTH(  J)#Jal/S03 _ 

STV(3)=STV(3)*((CALC( 17)-STV{2))*(CALC(17>-STV(2))) 


-STya3>sSTW<  3)*((CALC(  in)-STW(2))*CCALC.Cl8)-STVIC2XX)- 


STX(3)=STX(3)+((CALC( 16)-STX(2>)*(CALC(16)-STYC2) ) > 


STZ(3)  =  ST7(  3)«>((CALC(20)-ST7.(2))*(CALC{20)-STZ(2))) 

j»a-COffT-INUC - 

STZA(2)=STZA(2)/C0UNT 

- OO-LOS  llARKJtl^KOUNT - 

READ!  11>MARK/ 53  SMOOTH C.J),Jsl/r30) 

- S4'-2A^3)-=5TZA-(-3-H444«OOTH(1-)«ST.7A<2)>«<SHQOTHU>^TZJU2))) 

105  continue 

- C<JUHTVaKOUMT»'l - - - - - 


«RITE<  I01IT#76) 


lERRQR^lOX^iaHMEAN  FORCAST  ERROR/ 12X,flHVARIANCE////) 

___SIS<-t3  X.-  S-IVLCl  J  /  C  Qt?  N I J _ _ _ _ _ _ 

STVCl)  =  STV(l)/CniiNT 

_ iiRii£XiouT.7-n(;s-Tvcj)/j=Liii _ 

77  FORKAT(lHO/lOX/lAHMQWIHG  AVERAGE/ 18X/ 3(Fa.2/20X )////) 

eTu/1^-CTur^^/rn^tMT1 


STWCDsSTWd  WCQUNT 

_ laUE  C.I  OUT-/  zaj  ( S.T  w  c  J-i/  J-=  l.  ix _ _ _ 

7d  K0RHAT(1H0/10X/13HLEAST  SaUARES/19X/3(Fa.2/20X)////) 

_ ^TXCa)g3TXf  3>/CaUUT  l _ - 

STX(l)=STX(l)/COrNT 


iM*n  #10  A 


79  FQRHAT(lHO,5X,2aHSINr,LF  EXPONENTIAL  SH0nTHING/9X/3(F8.2/20X>////) 

- SWa^)-i3.TY<-3-)./CQUNT-t - 

STY<nsSTY(n/Cn(tNT 

- HR-ttE-d  niLT.-30  V(  ST  YC  J  )  - - 

30  FQRMATClH0»5X#28H0UUnLE  EXPONENTIAL  SMQ0THING»9X*3(r8.2/20X)////) 

\  >rn«» MT4 


STZCDsSTZCD/COUNT 

- NR-I-T£TI0UT/31->  (  STZCU)/JaI/3.) - 

31  FQRHAT(lH0/5X,2ftHTRlPLE  EXPONENTIAL  SMOOTHING  / 9X/ 3<Fa .2/20X)//// ) 

_ Si2AC3X-3-S.T2A(3WCQlJNTl - 

STZA(l)  s  ST7A(1)/C0UNT 


32  FORHATC 1H0/IX/37HA0APTIVE  SINGLE  EXPONENTIAL  SMOOTHING  /  4X/ 

_ UX£ia*2/20XJ////3 - 

HRlTE<lOUT/33) 

Ol-EORKAJ  ( IH 13 _ 

STOP 


5  0*14  P.H.  ASR  5.4  USING  73/345  XFORTN  COMPILER 

TIME _ 1 98  -SECS - 432-  CAROS.  A-TL— 1 32_C*P-*H . _ d.  f LA.( 

■  0  DATA  =  6174  TEMPORARIES  *  72  CODE  *  27900 
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